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STRESSES CAUSED BY AN IN-PLANE LINE LOAD USED TO 
ROTATE A TfflN- WALLED SPHERICAL SATELLITE 


By W. Jefferson Stroud and Nancy P. Sykes 
Langley Research Center 

SUMMARY 

The stress resultants and deformations which are brought about by a sinusoidal 
in-plane line load acting along a great circle of a spherical satellite of the Echo type are 
calculated and are presented in contour plots and in tabular form. The line load is 
applied in order to rotate the thin-walled spherical satellite for station-keeping purposes. 
The line load is induced by the interaction of a current-carrying wire with the geomag- 
netic field. Linear membrane theory and numerical methods are used in the analysis. 

INTRODUCTION 

One method for obtaining orbit position control of large, inflatable, spherical, pas- 
sive communications satellites of the Echo type is to utilize forces caused by solar photon 
flux on the satellite's surface. These forces are dependent upon the reflectance, absorp- 
tance, and emittance properties of the surface of the satellite. The satellite's surface 
could be divided into regions having different values of these properties. By adjusting 
the orientation of the satellite so that a desired surface faces the sun, the forces caused 
by solar radiation pressure could be modulated and used to control the orbit of, or 
"station keep," the satellite. To obtain the desired orientation the satellite could be 
rotated by torques resulting from the interaction of the earth's magnetic field with elec- 
tric current carried by wires mounted on the surface of the satellite. A discussion of 
this type of station-keeping system is contained in reference 1. 

In the present paper an analysis is made of the stresses in and deformations of a 
thin-walled spherical shell with a single current-carrying wire mounted on a great circle 
of the sphere and interacting with the geomagnetic field. The interaction causes both 
in-plane and radial loads on the sphere. Only the effects of the in-plane component of 
the applied line load and the resulting inertia load caused by angular acceleration are 
considered. Linear membrane theory and numerical methods are used in the analysis. 
Results are given as contour plots and as a table of stresses and displacements over the 
surface of the spherical shell. 



SYMBOLS 


The physical quantities used in this paper are given both in the U.S. Customary 
System of Units and in the International System of Units (SI) (ref. 2). Appendix A pre- 
sents factors relating these two systems of units. 

A = IB cos /3 

a^ constant associated with homogeneous solution of differential equation 

for Ng q ( see ec l- ( B7 )) 

B flux density of geomagnetic field (see eq. (1)) 

b direction of geomagnetic field (see figs. 1 and 2) 

D bending stiffness 

E Young's modulus 

dF differential force on wire 

dFj in-plane component of dF 

dF r out-of-plane (radial) component of dF 

I current in wire (see eq. (1)) 

I area moment of inertia of cross section of sphere (see eqs. (B13) and (B16)) 

Ip polar mass moment of inertia of sphere (see eq. (7)) 

K extensional stiffness 

dZ differential length along wire 

N n principal stress, for problem herein 

Ng, N^, in-plane stress resultants 
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N 6>,o> N e,i> — > N e,i> ••• 


functions of 9 which appear as coefficients in power 
series expansion for Ng (see eq. (B4)) 


^9<p O’ •••> ^94>,i’ ••• functions of 9 which appear as coefficients in power 

series expansion for (see eq. (C18)) 

P 0 ,P 0 ,P r loading per unit area in 9, <p, and radial directions, respectively 

r radius of sphere 

t wall thickness of sphere 


u,v,w displacements in 9, <p, and radial directions, respectively 

u 

U = T 

sin (p 

V shear load (see eq. (Bll)) 


v 0’ v l’"’’ v i’"' functions of 9 which appear as coefficients in power series expansion 
for v (see eq. (C19)) 

X,Y,Z axis system (see figs. 1 and 2) 


a 

P 


e 9> e <p>r9(p 


e,<P 


angular acceleration 

angle between plane of wire and earth's magnetic field (see figs. 1 and 2 ) 
in-plane strains 

circumferential and meridional coordinates of the spherical coordinate 
system (see fig. 2 ) 


A9,A(p increments in 9 and cp directions, respectively, used in the numerical 
calculation of N# 

i u Poisson's ratio 

| dummy variable of integration (see eq. (B14)) 
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p mass per unit area of sphere 

00 relaxation factor 


LOADING ON SPHERE 

The problem considered in this report is illustrated in figure 1, which shows a 
spherical shell with a current-carrying wire attached to the surface along a great circle. 
The arrows on the wire indicate the direction of current flow and the arrow at b 



Figure 1.- Thin-walled spherical satellite with current-carrying wire attached 
to surface along a great circle. 

indicates the direction of the geomagnetic field. The interaction of the current flowing 
in the wire with the geomagnetic field produces a force on the wire. The magnitude and 
direction of this force is given by the vector equation (see ref. 3) 

dF = T X B 61 (1) 

in which 

dF differential force, newtons 

T current in wire, amperes 
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B flux density of earth's magnetic field, webers per square meter 

dZ differential length along wire, meters 

The differential force dF can be divided into a component dFj which lies in the plane 
of the sphere and a component dF r which is directed normal to the surface of the 
sphere, or radially. The magnitude of the in-plane component dFj is given by the 
scalar relation 

dFi = IB dZ cos j8 sin cp (2) 

whereas the out-of-plane, or radial, component dF r is given by 

dF r = IB dZ sin j3 (3) 

An octant of the sphere of figure 1 is shown in greater detail in figure 2. Spherical 
coordinates <p and 6 are used to locate a point on the surface of the sphere. The 
meridional angle 4> is measured from the Z-axis, which passes through the upper pole. 
The circumferential angle 9 is measured from the X-Z plane, which is also the plane 
of the wire. Angle p is the angle between the plane of the wire and the direction of 
the geomagnetic field, axis b. Axis b lies in the X-Y plane so that angle p is mea- 
sured in the same way that angle 9 is measured. 


z 



Figure 2.- Octant of sphere with current-carrying wire showing differential loading 
and spherical coordinates. 
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If both the angle 0 and the direction of the current I flowing in the wire are as 
shown in figures 1 and 2, then the radial force dF r , which is independent of 4 b, is 
directed inward and the in-plane force dFj is directed in the 0 direction. 

In this report only the in-plane component of the loading is considered. The con- 
stant radial component is treated on page 338 of reference 4. Since linear theory is used 
in reference 4 and in the present report, the effects of the two loadings may be 
superimposed. 

The in-plane line load per unit length, dF^/dZ, introduced into the shell by the 
current-carrying wire may be written as 


dFi 

— — = A sin 4 > 


( 4 ) 


where, from equation (2), 


A = IB cos j3 


( 5 ) 


This applied sinusoidal line load (illustrated in fig. 3) gives a resultant torque 
shell of magnitude 

T = f Ar^sin^ d cp = Av^tt 
J 0 


T on the 
( 6 ) 


From Newton's second law, for rotational motion, 

T = I n a = — fl-rW 
p 3 


( 7 ) 



Figure 3.- Sinusoidal in-plane line load acting along a great circle of a sphere, 
causing a resultant applied torque. Top view. 
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Thus, the angular acceleration is 


a = 


3A 

8r 2 p 


( 8 ) 


The applied torque is equilibrated by an inertia loading which is distributed over the 
surface of the shell and which is given by 

Pg = -pra sin 0 (9) 

or, using the value of a just calculated, 

P e =-|£sin 0 (10) 


Thus, the in-plane line load (eq. (4)) and the surface inertia load (eq. (9)) have been 
expressed as a function of the strength of the geomagnetic field, the amount of current 
flowing in the wire, the orientation of the wire with respect to the geomagnetic field, the 
radius of the spherical shell, the mass density of the sphere, and the position on the 
sphere. 


SOLUTION OF MEMBRANE SHELL EQUATIONS 


Equilibrium Equations and Stress Boundary Conditions 


Because of the nature of the loading (all forces tangential to the surface) linear 
membrane shell theory is assumed to be adequate for investigating the stresses and 
deformations in the present problem. The coordinate system and a loaded shell element 
are shown in figures 2 and 4, respectively. The applied in-plane line load acts in the 
0-direction along the great circle 0=0 and 9 = tt. 



From reference 5 the equilibrium equations 
for a nonaxisymmetrically loaded spherical mem- 
brane are given by 

~W + sin $ + 2 N e<t> cos ^ + p 0 r sin <#> = 0 

(ID 

9Noj. 9N. 

- N 0 cos 0 + -^2. sin 0 + N^ cos (p 

+ P^r sin 0 = 0 (12) 

N0 + N^ = rP r (13) 
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These equations represent three equations in the three unknown stress resultants N$, 
N^, and which are to be solved in the region 0 ^ 0 = 77 , 0 S <p S it. For the prob- 
lem considered herein, P r = = 0; therefore, = -Ng. 

A complete set of boundary conditions can be written for 9 = 0,7 r and (p = 0,77. 
From symmetry conditions and from consideration of the applied loading given by equa- 
tion (4), boundary conditions on N 0 can be written along 9 = 0,77. The remaining con- 
ditions, those on (which is equal to -N$) and along < p = 0,77 and those on 

N along 0 = 0,77, are derived in appendix B. A summary of the stress boundary con- 
ditions is as follows: 


Along 6 = 0, 

N 0 = - | sin 4> 

(14) 


™ 6 <p = | Ac °s <f> 

(15) 

Along 9 = 77, 

N 0 = j sin (p 

(16) 


N 0 (p = | A cos (p 

(17) 

Along (p = 0, 

N^ = -N 0 = |Asin20 

(18) 


= | A cos 20 

(19) 

Along <p = 77, 

N0 = -N 0 = |Asin20 

(20) 


N 0 <p = - ^ A cos 20 

(21) 

In addition, symmetry conditions provide knowledge of certain stresses along 6 = ■=• 

77 ^ 

and (p = ± as follows: 


Along 6 = |, 

N s = 0 

(22) 

11 

■e- 

o 

3 

N M> = ° 

(23) 


Stresses 

Solution for the stress resultant Ng.- The governing partial differential equation 
obtained by eliminating N<£ and from the three equilibrium equations is 
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( 24 ) 


9 No 9 N fi / o o \ o 8 

^ + 5 sin <p cos <p + (4 cos^c p - 2 sin^JNa + sin^cp — = 0 

902 dfi \ > 9^,2 

With first and second derivatives replaced by two- and three-point central- 
difference approximations, respectively, equation (24) may be written in difference form 
as 


N e (e+A e,cj>) - 2N 0 (e,<#>) + N e (0-A0,4>) n 0 (0, <£+*<£) - N e (0,4>-A4>) 

(A 0)2 + 2A0 


5 sin 4> cos cp 


+ N^(0, 4>)\^ cos^<p - 2 sin^c^j 


N fl (0,4>+A<£) - 2N Je,<p) + N Je,<p-A<p) „ 

+ 6 L 6 sinrcp = 0 

(A^) 2 


(25) 


An iteration technique known as successive overrelaxation (ref. 6) is applied to 
equation (25) in the region 0 S 0 = ir/2, 0 S <p ^ tt. The boundary conditions on Ng are 
given in equations (14), (18), (20), and (22). Equation (25) written in successive over- 
relaxation form is 


N< n+I) (e,« = 


■00 


+ 4 cos 2 $ - 2 sin 2 <fi - 2 sin V 


/ ( 0+A 0, <p ) + N^ n+1) ( 0- A 0, (p ) 


J (A0)‘ 


(A0) 


(a <py 


N ( fl n) (0,<^>+A0) 


5 sin (p cos cp 
2 A <p 


sin 2 (ft 
( A ^) 2 


N ( fl n+1) (0,^-Arf>) 


5 sin 0 cos 4> si rr<t> 
2 A 0 (A 0) 2 


- (w - l)N< fl n) (0,<?>) 


(26) 


in which a> is the relaxation factor and the superscript n refers to the value of Ng 
obtained on the nth iteration. The purpose of the relaxation factor u> is to speed up 
convergence of the procedure. With co = 1.0 the method reduces to the better known 
but more slowly converging Gauss-Seidel iteration technique. The determination of an 
appropriate value for u> is discussed in reference 6. For the calculations performed 
in this study the value co = 1.95 was used. In the form given by equation (26), the com- 
putation is initiated at the point (0=A0, <p=A(p), which is one station to the right of the 
0 = 0 boundary, and proceeds in the 0 direction to the point ^0=^-A0, <p=A(pj, which is 
one station to the left of the 0 = ^ boundary. The coordinate <p is then increased by 
the increment A <p and the computation again proceeds from 0 = A0 to 0 = ^- A0. 
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This procedure is repeated until equation (26) has been applied to all points in the region 
0 < 8 < 1 ~, 0 < cf> < 7 r. In the present problem, a linear distribution in the 0 direction 

using the values of N q along the lines 0 = 0 and 0 = | was used as an initial approxi- 
mation to start the computation. The complete cycle is repeated until convergence is 
obtained. 

A contour plot of N 0 /A is shown in figure 5 and the values are given in table I. 

For more details concerning methods of solution and for estimates of the accuracy 
of all stress and displacement results, see appendix C. 


e 


0 


7T 

4 



7r 3n 

2 4 



Figure 5.- 


Contour plot of 


A 


in the hemispherical region 
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Solution for the stress resultant N g^.- Equation (11) can be used to determine 


N Q(p. Introducing into equation (11) the expression for P g given in equation (10) and 
dividing the equation by sin <fi yields 


8N 60 3A 1 8N 0 

90 8 ^ sin 0 00 


2N0<^ 


cos 
sin <p 


( 27 ) 


Equation (27) is integrated numerically in the <p direction for Ng<£ by using a Runge- 
Kutta method (ref. 7) for each value of 8. The derivative 9N g/ 90 is evaluated numeri- 
cally by using the Ng distribution stored in the computer. A contour plot of Ng^/A 
is shown in figure 6, and the values are given in table I. 
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Principal stress resultant N n .- By making use of the fact that = -N#, the 
principal stress resultant N n can be calculated from the equation 

N n = (n 0 2 + N e /) 1/2 (28) 

The values of the nondimensionalized principal stress resultant Nn/A are shown as a 
contour plot in figure 7 and are given in table I. 



Figure 7.- Contour plot of in the hemispherical region O = 0 = 7r,O = tf>=7r. 
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Displacements 


The displacements u, v, and w in the 9, and radial directions, respectively, 
are shown in figure 8. They are calculated by using the following strain- displacement 
relations (ref. 4, p. 96): 


3v 

— r + w = re , 
3 (p 4* 

(29) 

+ v cos 4> + w sin 4> = r eg sin <p 

(30) 

|H- sin 0 - u cos cf> + || = r Yg<j) sin <t> 

(31) 


The strains may be expressed in terms of the stresses by using 



The extreme right-hand sides of (32) and (33) were obtained by using = -N q. 



Figure 8.- Displacements u, v, and w in the 9, <P, and radial directions, 
respectively, of a point on the surface of a sphere. 
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Solution for the displacements v and w .- The governing partial differential 
equation for v is obtained from equations (29) to (34) and is given by 


3^v 


O^v • 2 j. 3v • . . 

— - sin 0 - — sm (p cos cp + v + — - 

3<£ 2 d( P 30 2 


-r sin ^4> 


±_m) 


Et 


^H-rsin ^ 2 (1 + M) 9 % 


3 <p 


Et 


90 


(35) 


By making use of equilibrium equations (11), (12), and (13) the right-hand side of (35) can 
be changed to give 


sin 2 0 - — sin 0 cos 0 + v + = - r - ^ + ^ sin 0 cos 0 (36) 

30 2 30 90 2 Et 


The boundary conditions for v, which may be obtained from symmetry and antisymmetry 
conditions, are as follows: 


Along 9 = 0, |, 
Along (p = 0, 


v = 0 


v = 0 


(37) 

(38) 


The displacement v is calculated in the same way that the stress resultant Ng 
is calculated — by the successive-overrelaxation method. The displacement w is cal- 
culated from equation (29). Contour plots of vEt/rA(l + p) and wEt/rA(l + p) are 
shown in figures 9 and 10, respectively, and the values are given in table I. 


Solution for the displacement u/sin cp.- Rather than solving for the displacement u 

in the 0 direction, it seems preferable to solve for u = -r . The displacement u is 
’ sin <p 

less useful than u because u would be zero at <p = 0 and <p = it even if a rigid-body 
rotation in the 0-direction were allowed. The quantity u/r is the angle of rotation in 
the 0-direction of a point on the shell as the shell deforms. At the pole cp = 0, the 
quantity u/r is the angle of rotation of a line element passing through the pole. The 
equation for u is obtained from equation (31) and is given by 


3u_ 

3 ( p 


3v 2r(l + jit) 
30 + Et 


^9(p 


sin 0 


sin ^0 


(39) 


The behavior of du/Bcp at the pole (p = 0 is investigated by expressing the 
numerator and denominator on the right-hand side of equation (39) in powers of (p . The 
distribution of u at (p = 0 is determined by expressing the terms in equation (30) in 
powers of <p. The resulting values of Bn/dcp and u at the pole are 
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(0 = 0 ) 


(40) 




H=-^^(sina + sin 3e) 

u = 3rA 8 ( E t + M) - cos 26 (0 = 0) (41) 

Because of the difficulty of accurately evaluating the right-hand side of equation (39) near 
(but not at) the pole, interpolation formulas are used in that region. See appendix C for 
details. 

With the above information it is possible to solve for u by using equation (39) and 
a Runge-Kutta method. A contour plot of uEt/rA(l + /i) sin 0 is shown in figure 11 and 
the values are given in table I. 



Figure 9.- Contour plot of in the hemispherical region 0 = 0 = 7r, 0 = cp ^ 7r. 

rA(l + m) 
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DISCUSSION OF RESULTS 


Contour plots of 


vEt 


wEt 


N 0 N g(/) N n 
A ’ A ’ A ’ rA(l + i±)’ rA(l + fi )’ 


and 


uEt 


are 


rA(l + p) sin <f> 

shown in figures 5, 6 , 7, 9, 10, and 11, respectively, and the values of these nondimen- 
sionalized stress resultants and displacements are given in table I. The inadequacy of 
linear membrane theory near the line loads is apparent in that the displacement w is 
discontinuous at the line loads 0 = 0 and 0 = tt. In the regions near the lines 0 = 0 
and 0 = 7T, bending occurs and/or nonlinear membrane action is present which, because 
of antisymmetry, would cause w to be zero. These linear membrane theory results 
are expected to be correct a short distance away from the line loads. 


Inasmuch as the station- keeping method referred to herein would be applied to very 
thin Echo-like satellites which carry no internal pressure, the compressive stresses 
caused by the line loads might bring about buckling of the surface. The magnitude of the 
load that would cause buckling can be approximated by making use of results of buckling 
tests on similar spheres presented in reference 8 . The uniform pressure differential 
that caused buckling in a 12.5-ft (3.81-m) diameter sphere of aluminum foil-Mylar 
laminate was about 14 percent of the classical buckling pressure differential of 


2 

Ap c i = (|) l/DK (42) 

where D is the bending stiffness and K is the extensional stiffness of the laminate. 
The stress resultant N n caused by this pressure differential is 

Ap rl r 

N n,cl = —f- (43) 

For the problem of a line load on a sphere considered in this report, the maximum value 
of the principal stress N n is A/2 and occurs at the points (0 = 0, <t> =Z ) and 

/ 7T\ 

( 0 = 77, (p - 2 )- If if is assumed that the surface of the sphere buckles at the same maxi- 
mum stress level under a nonuniform in-plane load as under uniform external pressure, 
then local buckling is initiated when 


A 0. 14(|f l/DKr 
iN n,max 2 ~ 2 

From equation ( 8 ), the angular acceleration a is related to the constant A by 


a = 


3A 

8r^p 


(44) 


( 45 ) 


The values of A and a which cause buckling are now calculated for the cases 
of a 12.5-ft (3.81-m) diameter sphere and a 135-ft (41.2-m) diameter sphere. For 
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each case the time t that is required to rotate the initially nonrotating sphere through 
an angle of n/2 radians is also calculated from the formula 

| = | o-t 2 (46) 

Note that the value of the angular acceleration a given by equation (45) causes buckling 
of the surface of the sphere; therefore, it is considered the upper limit for a. If this 
constant limiting value of a is used in equation (46), the value of the time t is the 
minimum time required to rotate the sphere through an angle of tt/2 radians. 

Both spheres are assumed to be fabricated of the aluminum foil-Mylar laminate 
discussed in reference 8. For this material, which is nominally identical to that used in 
Echo A- 12 (Echo II), the material properties are 

D = 23.8 X 10" 5 lbf-in. (2.69 X 10” 5 N-m) 

K = 3055 lbf/in. (53.50 X 10 4 N/m) 

p = 1.80 X 10"® slugs/in 2 (40.7 g/m^) 

The values for D and K are taken from reference 8. Results for the two spheres are: 

Case 1 Case 2 

Diameter = 12.5 ft (3.81 m) Diameter = 135 ft (41.2 m) 

r = 75 in. (1.90 m) r = 810 in. (20.6 m) 

A = 0.00637 lbf/in. (1.11 N/m) A = 0.00059 lbf/in. (0.103 N/m) 

a = 2.83 rad/sec^ a = 0.00225 rad/sec^ 

t = 1.05 sec t = 37.4 sec 

These results indicate angular- acceleration limitations caused by geometrical and 
material factors. For the second case, corresponding to the Echo A- 12 (Echo II) satel- 
lite, the minimum time required for rotation would not appear to impose any undue opera- 
tional restrictions. However, other limitations may be associated with concentrated 
masses such as battery packs, tracking beacons, and other instrumentation attached to 
the surface of the sphere. 


CONCLUDING REMARKS 

Stress resultants and displacements have been calculated for a sinusoidal in-plane 
line load acting along a great circle of a spherical passive communications satellite. 
Linear membrane theory is used and, within the framework of this theory, other loadings 
may be superimposed on the given one. An example is presented of a thin-walled spheri- 
cal satellite with a load applied which initiates local buckling at points on the surface of 
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the satellite. The results from this example indicate how the buckling strength of the 
structure may impose limitations on the maximum angular accelerations which can be 
applied to the satellite. Also, the results indicate that the order of magnitude of these 
limitations for a passive communications satellite of practical proportions is not 
prohibitive. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., August 7, 1967, 

124-08-06-02-23. 
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APPENDIX A 


CONVERSION OF U.S. CUSTOMARY UNITS TO SI UNITS 

The International System of Units (SI) was adopted by the Eleventh General 
Conference on Weights and Measures, Paris, October 1960, in Resolution No. 12 (ref. 2). 
Conversion factors for the units used herein are given in the following table: 


Physical quantity 

U.S. Customary 
Unit 

Conversion 

factor 

(*) 

SI Unit 

Force 

Length 

Mass 

lbf 

in. 

slugs 

4.448222 
0.0254 
14 593.90 

newtons (N) 
meters (m) 
grams (g) 


^Multiply value given in U.S. Customary Unit by conversion factor 
to obtain equivalent value in SI unit. 
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APPENDIX B 


THE DEVELOPMENT OF BOUNDARY CONDITIONS 


Equation (24) is a partial differential equation of the second order in both 0 and 
<p. The unknown in that equation is N g. It is also possible to manipulate equations (11), 
(12), and (13) to arrive at a similar equation with N g<£ as the unknown. In order to 
solve either of these two equations it is necessary to have four boundary conditions — 
two along 8 boundaries and two along <p boundaries. At present, consider equa- 
tion (24) with Ng as the unknown (repeated here for convenience): 


9 N 0 9N 0 / x 9^N e 

+ 5 sin (f) cos <fi + (4 cos*^ - 2 sin^$)Ng + sin^ = 0 (Bl) 

de 2 9 0 B<t> z 


The two boundary conditions which are given in equations (14) and (22) can be 
developed on the basis of symmetry conditions. These same two boundary conditions 
satisfy the requirement for two boundary conditions along 8 boundaries and are given 

by 


Ng = - | sin 4> 

(0 = 0) 

(B2) 

o 

11 

£ 

II 

(B3) 


Symmetry of Ng about the line <p - ^ can be used to satisfy one of the two boundary 
conditions along 0 boundaries. The boundary condition which is needed is the dis- 
tribution of N q at a pole 0 = 0 or 0 = tt. 

Let N# be expanded in a power series in 0 as 

Ng = Ng >0 + + N e>2 0 2 + . . . + N 0ji 0i + . . . (B4) 

in which the coefficients Ng^ are functions of 8. These coefficients can be obtained 
by substituting equation (B4) into equation (Bl) in which the trigonometric functions have 
also been expanded in powers of 4>. After coefficients of like powers of </> are grouped, 
equation (Bl) becomes 

n 8,o + 4N «,o + ( N e,i + 9N e,i)* + ( n «,2 + 16N «,2 - 6N e,oV 2 + ■•■=<> (B5) 

It is valid to set equal to zero the coefficients of each power of 0 in equation (B5). 

When this is done a set of ordinary differential equations is obtained, the first of 
which is 


N 0,O + 4N 0,O " 0 


(B6) 
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APPENDIX B 

The solution to equation (B6) is 

N 0 0 = a i sin 2 # + a 2 cos 20 (B7) 

The boundary condition given by equation (B2) is also expanded in a power series. That 
series together with equation (B3) provides the following boundary conditions for 

N 0>o = 0 (0 = 0) (B8) 

N 0 , O =0 (0 = |) (B9) 

The constant ^ is therefore zero, but a^ is undetermined. The distribution of 
at <p = 0 is thus known, but the amplitude of that distribution is not known. 

The two boundary conditions along c p boundaries which are required in order to 
obtain a solution to equation (Bl) are thus 

Ng = a^ sin 20 (<p = 0, ir) (BIO) 

Suppose one arbitrarily assumes that aj is zero, uses that as the boundary condi- 
tion at the upper and lower pole, and solves for the stress distribution Ng using succes- 
sive overrelaxation in the region 0^0 = ^, 0 S <p S n. Three examples of the results, 

Ci 

which depend upon the grid size in the (p direction, are shown in figure 12. In each of 
the three cases the results exhibit wide oscillations very near the pole <£ = 0, but four 


A CD 



o H 27T 37T 47T 57T 67r 77T 87T 97T 10_7T 

80 80 80 80 80 80 80 80 80 80 

0 ) 


Nn 

Figure 12.- The nondimensionalized stress near 0) = 0 along the line 6 = J 

A an 4 

for three values of A<t>. For all three cases = 0. 
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or five stations away from the pole the solution converges to the correct value of N# 
regardless of the choice of A 0. This phenomenon is demonstrated again in figure 13, 
which shows the stress near the pole for three choices of aj. These results suggest 
that the smoothness of the Ng distribution near 0 = 0 can be used as a criterion for 
determining the correct value for aj. 

The value of a j can be determined to three places by using graphs similar to 
those shown in figures 12 and 13. More accuracy can be obtained by using several values 
of a^ and calculating a table of differences for each value of aj used. The dif- 

ferences on Ng can be taken in the ^-direction along the line 8 = ^ in the neighbor- 
hood of the upper pole. By using smoothness in the differences as the criterion, the 
constant a^A was determined to be -0.37503 ± 0.00004. The uncertainty (±0.00004) 
arises because the digit in the last place is affected by the mesh size used to converge 
on the constant. That this value of a j is essentially correct can be ascertained by 
another method. 

A top view of the loaded sphere divided along the line of action of the applied line 
load is shown in figure 14. Half the line load is assumed to be acting on the left hemi- 
sphere, and half on the right. The in-plane inertia load is represented by the arrows 



)^0 

Figure 13.- The nondimensionalized stress near 0 = 0 along the line 0 = - for 

a l 

three values of — 

A 
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Figure 14.- Top view of loaded sphere divided along the line of action of the applied 

line load. 

around the circumference. By considering the vertical* equilibrium of the hemisphere 
on the right, for example, it is possible to calculate the distribution of the shear stress 
along the cut which passes through the upper and lower poles. This knowledge provides 
the required boundary condition. 

The vertical* resultant of the inertia load can be considered to be a shear load V 
given by 

V = ^ ^ (p sin 9 d0 dcp (Bll) 

= | Aru (B12) 

O 

The shear-stress distribution along the cut can be related to this shear load by using the 
familiar equation for the transverse shear-stress distribution in engineering beam theory: 

br = 2N 6(j) = (B13) 

where b = 2t, t is the shear stress, and the cross section of the sphere is regarded 
as the cross section of an elementary beam. 


^Vertical here means the vertical direction in figure 14, which is the x-direction. 
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Neif 


The right-hand hemisphere of figure 14 
is shown in a view from center out in fig- 
ure 15. The large arrow indicates the direc- 
tion of the resultant shear force V. By- 
using this figure it is seen that if the shear 
stress is to be calculated at some 

point 0, then Q is the first moment of the 
material to the right of 0, the first moment 
being calculated with respect to the line 
0 = 0. Thus 

* 7T / 2 


Q=2 r 


0 


tr^sin £ d£ 


= 2tr^cos 0 


(B14) 

(B15) 


Figure 15.- Hemisphere showing resultant shear 
force V caused by in-plane inertia load Pq 
and the equilibrating shear stress N^. 


The quantity I, the area moment of inertia 
of the cross section of the sphere about the 
line 0 = 0, is given by 

(B16) 


I = 7rr 3 t 


By substituting equations (B12), (B15), and (B16) into (B13), the shear stress distribution 
along the line 0 = 0 is found to be 

3A 
8 


N 


00 


COS 0 


(B17) 


A shell element at the pole of the sphere, together with the positive directions for 
the stress resultants, is shown in figure 16. The values of the stress resultants along 


the line 0 = 0 are known: 


N 0 0 = ^ cos 0, 


Ng = - ■=- sin <j), and, in addition, 


N <p = 


-N 0 


everywhere (from eq. (13)). Therefore, the stress resultants for an element in the 
orientation shown in figure 16 are known. By using the transformation equations for 

plane stress it is possible to obtain the stress 
at this same point for other orientations - that 
is, for other values of 0. The stress resultants 
as a function of 0 at the point 0 = 0 are 
given by 



No = - — A sin 20 
U 8 

(B18) 

Na = ^ A sin 20 
' 8 

(B19) 

N 00 = ! A cos 20 

(B20) 
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A summary of the stress boundary conditions is 


N q = - ^ A sin 2 6 

(0 = 0 ) 

(B21) 

N 0 <£ = | A cos 20 

(0 = 0 ) 

(B22) 

N Qcj) = 0 

Nog 

II 

(B23) 

N 0 = - | sin 

(e = o) 

(B24) 

N# = 0 

II 

coN 

(B25) 

®e<p = |a cos 0 

(0 = 0 , 3T) 

(B26) 
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ACCURACY OF SOLUTION 
Stress Resultant Ng 

The distribution for Ng was calculated by using a five-point central-difference 
equation modified to make use of a successive overrelaxation technique. That difference 
equation (26) is repeated here: 


N ( 0 n+1) (^) = 


-0J 


N (n) (e+Ae,<M + 

2 + 4 cos 2 0 . 2 sin 2 ^, . 2 sln^\L * 


1 




J(a ey 




+ N ( e n) (0,0+A0) 


o 

5 sin (f) cos <p + sin 0 
2A <p (A 0) 2 


+ N ( e n+1) (0,0-A0) 


5 sin 0 cos 0 sin^cfr 


2A0 


(A0)‘ 


(w - DN ( e n) (0,0) 


(Cl) 


The boundary conditions are 


Ng = - £ sin 0 


Ng = 0 

Ng = - M sin 20 


(0 = 0) (C2) 

(C3) 

(0 = 0, tt) (C4) 


fl - 77 
e ~2 


The overrelaxation factor w was approximated as 1.95. The mesh size used in the 
final computation was A0 = t~~- and A <p = . 

The degree of accuracy that could be expected from the numerically obtained Ng 
distribution was ascertained by varying the mesh size in both the 0 and <f> directions 
and noting the digits which changed with the mesh size. Except near the pole ^0 < 0 < ^ 
the changes in Ng occurred in the sixth place. The results for Ng presented in 
table I are believed to be correct to the number of digits retained. 


Shear-Stress Resultant Ng^ 

The distribution of the shear-stress resultant Ng^ was calculated from the 
equilibrium equation (11), rewritten here in a slightly different form: 
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dcp 8 sin (p 8 0 


2N 


cos 4> 
W sin cp 


The following boundary conditions are known: 

N ecp = ^Y cos 26 


N 00 = 0 


(<P = 0) 

(♦■S 


(C5) 


(C6) 

(C7) 


The boundary condition (C7) was used to initiate the computation. Equation (C5) was 
integrated in the <£ direction by using the following Runge-Kutta method (ref. 7): 

+ 2kg + 2kg + k^ 

^n+l _ + g 



(C8) 


where 


ki = hf (x n , y n ) 

^2 = kf ^ x n + 2’ ^n + "2") 
k 3 = hf ^ n + |, yn +_ ^j 


k 4 = hf (x n +h, y n +k 3 ) 


In the notation of the present problem, 


y = N e0 
f 

dcp 

X = <p 

h = 2A (p 


where A cp is the mesh size used in the N e calculation. The integration progressed 
from <P =x toward the pole. 

A 

Because it appeared that the error in this integration method was beginning to build 

7T 

up as the integration neared the upper pole, the values of N q ^ in the region 0 < 6 = 

— S <p S — were calculated by using equation (12), which, for the present problem, 

20 20 

reduces to 

8N flA 9N fl 

— = 2No cos cp + - sin <p (C9) 
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The same Runge-Kutta method which was used to integrate equation (C5) in the 
0-direction was also used to integrate equation (C9) in the ^-direction. 

By varying the mesh size it was determined that the values of N qq as presented 
in table I are in error by no more than 1 unit in the last place retained. 

The Displacement v 

The displacement v in the 0-direction was calculated by the same method as that 
used to calculate Ng, a five-point central-difference equation utilizing successive over- 
relaxation. The difference equation is 



4r ^ E * — Ng(g,0)sin 0 cos 0> - ( w - 1 )v^(0,0) 


(CIO) 


v = 0 (0=0) (Cll) 

v = 0 (0 = |) (Cl 2) 

v = 0 (0 = 0) (C13) 

v = 0 (0 = |) (C14) 

Because the right-hand side of the partial differential equation for v, equation (35), 
could be modified in such a way that Ng is the only unknown (eq. (36)), the mesh size 
could be made as fine for v as it was for Ng. For this reason, good accuracy was 
obtained in the calculation of the displacement v. Thus, the values of v given in 
table I are correct to one unit in the last digit retained. 


The Displacement w 

The displacement w normal to the surface of the sphere was calculated from the 
strain-displacement relation (29). That relation, written in terms of differences and 
with the stress resultant Ng replacing the strain, is 


w * - r ( 1 - H) 
Et 


_ v(0,0+A0) - v(0,0-A0) 


(Cl 5) 
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The results of varying the mesh size indicate that w is correct to four places to the 
right of the decimal point. That number of digits is retained in table I. 


The Displacement u 


As indicated in the body of the report, rather than calculating the displacement u 
in the 0-direction, the quantity u = ^ was calculated. The equation used for the 

calculation is 


9u_ 

9 ( p 


9v 2r(l + u) 
90 Et 


N 0( £ sin cj> 


sin 2 0 


(C16) 


The behavior of this expression at <p = 0 is determined by expressing the numera- 
tor and denominator on the right-hand side of equation (C16) in powers of <p. In order 
to do this it is first necessary to have the series expansions for Ng, N^, and v. 

Let 

N0 = N 0>o + N e>1 0 + N^ 2 + . . . (C17) 

N 00 = N 00,O + N 0</>,10 + n 00,2^ 2 + • • • (C18) 

v = Vq + Vj0 + v 2 0 2 + . . . (C19) 


The second term in the Ng series is obtained in the same way that Ng g 

was obtained in appendix B - that is, by substituting the series expansion for Ng into 
the governing partial differential equation, equation (Bl), by setting equal to zero the 
coefficients of each power of 0, and by utilizing boundary conditions which have also 
been expanded in powers of cf>. The first two coefficients are 

N0, O = - ^ sin 20 (C20) 

n M = ' f cos 30 (C21) 


The first two terms in the series expansion for Ng^ are obtained in a similar 
manner. All terms in equation (11) are expanded in powers of cf> and the equation is 
assumed to hold for each power of 0. By using the results given in equations (C20) and 
(C21), the quantities Ng^g and N g^i are given by 

n«*,o - x cos 28 <C22) 


N 00,l 


- ^ sin 30 

Ci 


(C23) 
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In the same way, by using equation (36), 


o 

II 

o 

> 

(C24) 

V 1 = 3Ar ‘ 1 Et + Sta 20 

(C25) 

v 2 = Ar 4Et ~( cos 36 ‘ cos 

(C26) 


After all pertinent results are substituted into equation (Cl 6) written as 


dV l 2r(l + u) XT 
9u = L JjT^ Et N e0,oj 0 + 

d 4> 6 


. 2r(l + m) 






it is seen that there is no singularity at the pole and 


9“ = - £A(1 t JL) (sin 0 + sin 30) 

9 <P 0=0 4Et 


By expanding u in powers of 0 as was done with N$, N^, and v, the distribution 
of u at 0 = 0 is determined as 

u(0,O) = + M)cos 20 + c (C29 

oEt 


where c is a constant which is affected by a rigid-body rotation. In the present analy- 
sis c is set equal to zero. 

Because of the small differences that occur in the numerator of equation (C16) for 
values of 0 near 0 = 0, the value of 9u/9 0 adjacent to the pole was calculated with 
an eight -point interpolation formula utilizing the value of 9u/9 0 calculated at the pole 
and seven other values calculated at points at least two stations away from the pole. 
Equation (C28) was used to obtain the value of 8u/9 0 at the pole, and (C16) was used 
to calculate 9u/9 0 at stations not adjacent to the pole. 

With the above information it is possible to solve for u by using equation (C16) 

uEt 

and a Runge-Kutta method. A contour plot of is shown in figure 11 and 

rA(l + p)sin 0 

the values are given in table I. 
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TABLE I.- NUMERICAL VALUES OF THE STRESS RESULTANTS AND DISPLACEMENTS 

= radians] 



N 0 


uEt 

vEt 

wEt 

N n 


n 

N 0 

Nflj, 

uEt 

vEt 

wEt 

Nn 


A 

A 

rA(l + ii)sin <p 

rA(l + f±) 

rA(l + jjl) 

A 


A 

A 

rA(l + fi) sin (p 

rA(l + fi) 

rA(l + fi) 

A 




9-0 








3?r 

0 = 2 ^- radians 



0 

0.00000 

0.3750 

0.375 

0.00000 

0.0000 

0.3750 


0 

-0.30341 

0.2204 

0.220 

0.00000 

0.0000 

0.3750 

1 

-.07822 

.3704 

.375 

.00000 

.0782 

.3796 


1 

-.30951 

.1517 

.169 

.04280 

.0697 

.3447 

2 

-.15451 

.3566 

.375 

.00000 

.1545 

.3889 


2 

-.30702 

.0981 

.125 

.07498 

.1377 

.3223 

3 

-.22700 

.3341 

.375 

.00000 

.2270 

.4040 


3 

-.30005 

.0587 

.087 

.09595 

.2022 

.3057 

4 

-.29389 

.3034 

.375 

.00000 

.2939 

.4224 


4 

-.29127 

.0309 

.056 

.10587 

.2618 

.2929 

5 

-.35355 

.2652 

.375 

.00000 

.3536 

.4420 


5 

-.28232 

.0125 

.031 

.10552 

.3150 

.2826 

6 

-.40451 

.2204 

.375 

.00000 

.4045 

.4607 


6 

-.27421 

.0016 

.011 

.09606 

.3604 

.2742 

7 

-.44550 

.1702 

.375 

.00000 

.4455 

.4769 


7 

-.26753 

-.0036 

-.004 

.07898 

.3969 

.2676 

8 

-.47553 

.1159 

.375 

.00000 

.4755 

.4894 


8 

-.26259 

-.0046 

-.014 

.05601 

.4236 

.2626 

9 

-.49384 

.0587 

.375 

.00000 

.4938 

.4973 


9 

-.25958 

-.0030 

-.021 

.02902 

.4400 

.2596 

10 

-.50000 

.0000 

.375 

.00000 

.5000 

.5000 


10 

-.25856 

.0000 

-.023 

.00000 

.4454 

.2586 




6 = ~ radians 







0 = |^ radians 



0 

-0.11589 

0.3567 

0.357 

0.00000 

0.0000 

0.3750 


0 

-0.35667 

0.1159 

0.116 

0.00000 

0.0000 

0.3750 

1 

-.18050 

.3196 

.334 

.01734 

.0773 

.3670 


1 

-.33091 

.0535 

.061 

.04905 

.0633 

.3352 

2 

-.23465 

.2786 

.313 

.03199 

.1526 

.3643 


2 

-.30500 

.0110 

.016 

.08415 

.1250 

.3052 

3 

-.27904 

.2377 

.294 

.04280 

.2242 

.3666 


3 

-.28132 

-.0159 

-.021 

.10587 

.1836 

.2818 

4 

-.31467 

.1978, 

.277 

.04905 

.2903 

.3717 


4 

-.26088 

-.0310 

-.051 

.11525 

.2377 

.2627 

5 

-.34261 

.1599 

.262 

.05043 

.3492 

.3781 


5 

-.24394 

-.0373 

-.075 

.11363 

.2860 

.2468 

6 

-.36390 

.1242 

.250 

.04708 

.3995 

.3845 


6 

-.23048 

-.0369 

-.094 

.10258 

.3272 

.2334 

7 

-.37946 

.0908 

.240 

.03946 

.4400 

.3902 


7 

-.22029 

-.0316 

-.109 

.08382 

.3604 

.2225 

8 

-.39003 

.0594 

.234 

.02836 

.4696 

.3945 


8 

-.21318 

-.0229 

-.119 

.05919 

.3847 

.2144 

9 

-.39615 

.0293 

.230 

.01482 

.4877 

.3972 


9 

-.20899 

-.0120 

-.124 

.03059 

.3995 

.2093 

10 

-.39816 

.0000 

.228 

.00000 

.4938 

.3982 


10 

-.20760 

.0000 

-.126 

.00000 

.4044 

.2076 




9 ~ — radians 
20 







0 = radians 



0 

-0.22044 

0.3034 

0.303 

0.00000 

0.0000 

0.3750 


0 

-0.37500 

0.0000 

0.000 

0.00000 

0.0000 

0.3750 

1 

-.25896 

.2433 

.263 

.03199 

.0744 

.3554 


1 

-.32419 

-.0432 

-.050 

.05044 

.0553 

.3271 

2 

-.28476 

.1895 

.226 

.05741 

.1470 

.3420 


2 

-.28273 

-.0675 

-.090 

.08503 

.1092 

.2907 

3 

-.30121 

.1441 

.195 

.07498 

.2159 

.3339 


3 

-.24968 

-.0786 

-.123 

.10552 

.1605 

.2618 

4 

-.31109 

.1068 

.168 

.08415 

.2795 

.3289 


4 

-.22374 

-.0803 

-.149 

.11363 

.2078 

.2377 

5 

-.31655 

.0770 

.145 

.08503 

.3362 

.3258 


5 

-.20369 

-.0753 

-.171 

.11110 

.2500 

.2172 

6 

-.31922 

.0537 

.127 

.07826 

.3847 

.3237 


6 

-.18854 

-.0654 

-.187 

.09966 

.2860 

.1996 

7 

-.32027 

.0356 

.113 

.06487 

.4237 

.3222 


7 

-.17750 

-.0519 

-.200 

.08106 

.3150 

.1849 

8 

-.32050 

.0215 

.103 

.04627 

.4522 

.3212 


8 

-.16999 

-.0359 

-.209 

.05706 

.3362 

.1737 

9 

-.32044 

.0101 

.097 

.02406 

.4696 

.3206 


9 

-.16564 

-.0184 

-.214 

.02944 

.3491 

.1667 

10 

-.32039 

.0000 

.095 

.00000 

.4755 

.3204 


10 

-.16421 

.0000 

-.216 

.00000 

.3535 

.1642 
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TABLE I.- NUMERICAL VALUES OF THE STRESS RESULTANTS AND DISPLACEMENTS - Concluded 

[* = i radians] 


n 

N 0 

N 6><j> 

A 

uEt 

vEt 

wEt 

N n 


n 

n e 

N M 

uEt 

vEt 

wEt 

Nn 

A 

rA(l + ju)sin 0 

rA(l + (i ) 

rA{l + i±) 

A 


A 

A 

rA(l + ja)sin <£ 

> 

1 

rA(l + y) 

A 




677 

6 - — radians 
20 







Q<rf 

9 = — radians 
20 



0 

-0.35667 

-0.1159 

-0.116 

0.00000 

0.0000 

0.3750 


0 

-0.11589 

-0.3567 

-0.357 

0.00000 

0.0000 

0.3750 

1 

-.29210 

-.1317 

-.156 

.04708 

.0459 

.3204 


1 

-.08712 

-.2972 

-.363 

.01482 

.0122 

.3097 

2 

-.24429 

-.1344 

-.188 

.07826 

.0908 

.2788 


2 

-.06828 

-.2500 

-.372 

.02406 

.0241 

.2591 

3 

-.20879 

-.1292 

-.214 

.09606 

.1334 

.2455 


3 

-.05552 

-.2114 

-.381 

.02902 

.0355 

.2186 

4 

-.18237 

-.1184 

-.235 

.10258 

.1727 

.2175 


4 

-.04667 

-.1776 

-.389 

.03059 

.0460 

.1837 

5 

-.16278 

-.1038 

-.253 

.09966 

.2078 

.1930 


5 

-.04045 

-.1464 

-.397 

.02943 

.0553 

.1519 

6 

-.14842 

-.0862 

-.266 

.08898 

.2377 

.1716 


6 

-.03607 

-.1165 

-.404 

.02609 

.0633 

.1220 

7 

-.13818 

-.0664 

-.277 

.07212 

.2618 

.1533 


7 

-.03304 

-.0872 

-.409 

.02104 

.0697 

.0933 

8 

-.13134 

-.0451 

-.284 

.05064 

.2795 

.1389 


8 

-.03106 

-.0581 

-.413 

.01472 

.0744 

.0659 

9 

-.12741 

-.0228 

-.288 

.02609 

.2902 

.1294 


9 

-.02994 

-.0291 

-.416 

.00757 

.0772 

.0417 

10 

-.12613 

.0000 

-.290 

.00000 

.2938 

.1261 


10 

-.02957 

.0000 

-.417 

.00000 

.0782 

.0296 




7 77 

0 = — radians 
20 







0 = radians 

20 



0 

-0.30341 

-0.2204 

-0.220 

0.00000 

0.0000 

0.3750 


0 

0.00000 

-0.3750 

-0.375 

0.00000 

0.0000 

0.3750 

1 

-.23862' 

-.2061 

-.247 

.03946 

.0355 

.3153 


1 

.00000 

-.3090 

-.379 

.00000 

.0000 

.3090 

2 

-.19347 

-.1878 

-.270 

.06487 

.0701 

.2696 


2 

.00000 

-.2578 

-.385 

.00000 

.0000 

.2578 

3 

-.16144 

-.1679 

-.289 

.07898 

.1030 

.2329 


3 

.00000 

-.2168 

-.392 

.00000 

.0000 

.2168 

4 

-.13843 

-.1467 

-.305 

.08382 

.1334 

.2017 


4 

.00000 

-.1814 

-.400 

.00000 

.0000 

.1814 

5 

-.12181 

-.1244 

-.319 

.08106 

.1605 

.1741 


5 

.00000 

-.1491 

-.407 

.00000 

.0000 

.1491 

6 

-.10987 

-.1010 

-.330 

.07212 

.1836 

.1492 


6 

.00000 

-.1184 

-.413 

.00000 

.0000 

.1184 

7 

-.10149 

-.0766 

-.338 

.05831 

.2022 

.1272 


7 

.00000 

-.0885 

-.418 

.00000 

.0000 

.0885 

8 

-.09595 

-.0515 

-.344 

.04088 

.2158 

.1089 


8 

.00000 

-.0589 

-.422 

.00000 

.0000 

.0589 

9 

-.09278 

-.0259 

-.347 

.02104 

.2242 

.0963 


9 

.00000 

-.0295 

-.424 

.00000 

.0000 

.0295 

10 

-.09176 

.0000 

-.349 

.00000 

.2270 ; 

.0918 


10 

.00000 

.0000 

-.425 

.00000 

.0000 

.0000 


0 ~ ^ radians 


0 

-0.22044 

-0.3034 

-0.303 

0.00000 

0.0000 

0.3750 

1 

-.16853 

-.2623 

-.318 

.02836 

.0241 

.3118 

2 

-.13375 

-.2265 

-.333 

.04627 

.0477 

.2631 

3 

-.10979 

-.1952 

-.346 

.05601 

.0701 

.2240 

4 

-.09297 

-.1662 

-.357 

.05919 

.0908 

.1905 

5 

-.08102 

-.1383 

-.367 

.05706 

.1092 

.1603 

6 

-.07255 

-.1109 

-.376 

.05064 

.1250 

.1325 

7 

-.06667 

-.0834 

-.382 

.04088 

.1376 

.1068 

8 

-.06280 

-.0558 

-.387 

.02863 

.1469 

.0840 

9 

-.06060 

-.0279 

-.390 

.01472 

.1526 

.0667 

10 

-.05989 

.0000 

-.391 

.00000 

.1545 

.0599 


NASA-Langley, 1968 31 L-5088 


35 


FIRST CLASS MAIL 


National Aeronautics and Space Administration 

WASHINGTON, D. C. 
OFFICIAL BUSINESS 


POSTAGE AND FEES PAID 
NATIONAL AERONAUTICS AND 
SPACE ADMINISTRATION 


05U 001 56 51 30 S 68074 00903 

AIR FORCE WEAPONS L AB0R A TORY/ AF WL / 
KIRTLAND AIR FORCE BASE, NEW MEXICO 87117 


ATI' MISS MADELINE F. CANOVA, CHIEF TtCHNI 
i. I BRAKY /Wi It/ 


POSTMASTER: 


If Undeliverable (Section 158 
Postal Manual) Do Not Return 


"The aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof.” 

— National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and technical information considered 
important, complete, and a lasting contribution to existing knowledge. 

TECHNICAL NOTES: Information less broad in scope but nevertheless of 
importance as a contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: Information receiving limited distribu- 
tion because of preliminary data, security classification, or other reasons. 

CONTRACTOR REPORTS: Scientific and technical information generated 

under a NASA contract or grant and considered an important contribution to 
existing knowledge. 

TECHNICAL TRANSLATIONS: Information published in a foreign 
language considered to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information derived from or of value to NASA 

activities. Publications include conference proceedings, monographs, data 
compilations, handbooks, sourcebooks, and special bibliographies. 

TECHNOLOGY UTILIZATION PUBLICATIONS^ Information on tech- 
nology used by NASA that may be of particular interest in commercial and other 
non-aerospace applications. Publications include Tech Briefs, Technology 
Utilization Reports and Notes, and Technology Surveys. 


Details on the availability of these publications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION DIVISION 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 



